Dynamics of Spin Relaxation near the Edge of Two-Dimensional Electron Gas 



O 

o 

(N 

m 



Yuriy V. Pershin 

Center for Quantum Device Technology, 

Department of Physics, Clarkson University, 

Potsdam, NY 13699-5721, USA 

We report calculations of spin relaxation dynamics of two-dimensional electron gas with spin- 
orbit interaction at the edge region. It is found that the relaxation of spin polarization near the 
edge is more slow than relaxation in the bulk. That results finally in the spin accumulation at 
the edge. Time dependence of spin polarization density is calculated analytically and numerically. 
The mechanism of slower spin relaxation near the edge is related to electrons reflections from the 
boundary and the lack of the translation symmetry. These reflections partially compensate electron 
spin precession generated by spin-orbit interaction, consequently making the spin polarization near 
the edge long living. This effect is accompanied by spin polarization oscillations and spin polarization 
transfer from the perpendicular to in-plane component. 
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PACS numbers: 71.70.Ej, 72.15.Lh, 76.60.Es 



The prospects for creating a semiconductor-based spin- 
tronic device [E Hi IM Hi IM ly] have generated an emphasis 
on the studies of properties of electron sp in p olarization 
in semiconductor nanostructures [3) IE 0> HE Hll HE HE 
HJ, llE HE H3 ■ Great interest has been expressed in dy- 
namics of electron spin polarization [3, IE IE HE HH HE| ■ 
For instance, spin relaxation dynamics has been stud- 
ied in two-dimensional electron gas (2DEG) 0, IE| , two- 
dimensional channels (EHEL open Sinai billiards (2DEG 
with a lattice of antidots) |ll| , and ballistic quantum dots 
[l2| . It was found that the sarnple geometry [Il| as well 
as specific initial conditions 0, lal could have a significant 
effect on electron spin relaxation. 

D'yakonov-Perel' (DP) spin relaxation mechanism 1^ 
is the leading spin relaxation mechanism in many impor- 
tant experimental situations. In the framework of DP 
theory, initially homogeneous electron spin polarization 
exponentially relaxes to zero (or to some finite equilib- 
rium value) with time. However, DP theory was formu- 
lated for the bulk of a sample. Considering electron spin 
relaxation near the edge of 2DEG, one would expect the 
same relaxation scenario. This expectation, however, is 
not correct. We demonstrate in this Letter that the spin 
relaxation dynamics near the edge is rather unusual and 
can not be described by a simple exponential law, as 
follows from the DP theory. Wc observe a longer spin re- 
laxation time near the edge, spin polarization oscillations 
and spin polarization transfer from the perpendicular (to 
2DEG) to in-plane component. 

Let us consider a two-dimensional electron gas with the 
Rashba spin-orbit interaction HE] , which couples electron 
space and spin degrees of freedom: 



Hn = ah ^ {a-^Py - cFyPx) , 



(1) 



where a is a constant, a is the Pauli matrix vector corre- 
sponding to the electron spin, and p is the momentum of 
the electron confined in two-dimensional geometry. From 
the point of view of electron spin, the effect of spin-orbit 
interaction can be regarded as an effective magnetic field 



acting on electron spin. Momentum scatterings reorient 
the direction of this field, thus leading to average spin 
relaxation (DP relaxation). Intuitively, electron reflec- 
tions from the 2DEG edges should increase spin relax- 
ation time, since electron spin rotations will be partially 
compensated. In order to study this phenomena, we use 
both analytical and numerical approaches. 

Let us consider a two-dimensional electron gas occupy- 
ing half-space a; > (Fig. 1) with electron spins homoge- 
neously polarized in z direction, perpendicular to 2DEG, 
in the initial moment of time t — Q. From the experimen- 
tal point of view, the initial electron spin polarization 
as well as detection of spin polarization can be realized 
using the optical technique [2E| or any other appropriate 
method. It is convenient to describe the spin state of i-th 
electron via the spin polarization vector P^ = Tr( pia), 
where pi is the single-particle density matrix |2]| . In 
the framework of drift-diffusion approximation, the ini- 




FIG. 1: Schematic of the system: 2DEG occupies half-plane 
a; > 0. The influence of the boundary located at x = on the 
electron space motion is taken into account by introducing 
the mirror source S' . 



tial dynamics of the spin polarization density P, which is 
defined via P(Pr — ^ P,j (here (Pr — dxdy is a space 

volume), can be calculated using the following expression 



P{x,y,t) = 



G{x 



x',y' 



y\t)P\, 



X;y)-{x',y 



)dy'dx' , 

(2) 
where G{x — x' ,y — y', t) is the diffusion Green function 
(solution of a diffusion equation with a point source), and 
P' {x,y),{x',y') represents a contribution of the initial spin 
polarization density at point {x' ,y') to P{x^y,t). The 
structure of Eq. Q can be easily understood. Electron 
spin polarization density in a space volume with coor- 
dinates {x, y) at a selected moment of time t is given 
by a sum of spin polarization vectors of all electrons 
located in this volume. The diffusion Green function 
G{x — x' ,y — y' ,t) gives probability for electrons to dif- 
fuse from the point {x\y') to (a;,j/), while P' (x,y).{x',y') 
describes the spin polarization of these electrons. 

We note that Eq. Q governs only the initial spin 
relaxation dynamics. The main approximation behind 
Eq. |(5J is the assumption that different spin rotations 
commute with each other, and spin precession angle ip is 
proportional to the distance between {x',y') and {x,y). 
This assumption is justified when the spin precession an- 
gle per mean free path is small and the time is short. 
Moreover, it is assumed that evolution of electron spin 
degree of freedom is superimposed on the space motion 
of the charge carriers. In other word, the influence of 
the spin-orbit interaction on electron space motion is ne- 
glected. If a is the unit vector along the precession axis, 
then ,]jj 
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FIG. 2: (Color online) Initial evolution of spin polarization 
near the edge calculated within the drift-diffusion approxima- 
tion. Relaxation of electron spins occurs slowly near the edge. 
All graphs presented in this Letter were obtained using the 
same parameter yLp = 0.1. 



where Pj_ = P — a(aP) is the component of the spin 
polarization perpendicular to the precession axis, Lp — 
Tjr, 77 is the spin precession angle per unit length, r — 
{x — x',y — y'), r = |r|, a = ^ x r/r, and z is the unit 
vector in z direction. We emphasize that the Rashba 
spin-orbit interaction is the origin of spin polarization 
rotations described by Eq. (Q. The coupling constant a 
of the Rashba interaction enters into our model through 
the parameter 77 (for details see Ref. 10]). 

When 2DEG occupies only a half-space, the Green 
function G{x — x', y — y', t) that appears in Eq. ^ can be 
obtained using mirror-image approach. Within this ap- 
proach, the Green function G{x — x' ,y — y' ,t) is written 
as a superposition of two point-source solutions 



G = Go{x -x',y- y\ t) + Go(x -x',y + y' , t), (4) 



Lp is the mean free path, and Tp is the momentum relax- 
ation time. While the first term in the right hand side 
in Eq. Q describes the real source, the second term is 
a so-called mirror source (Fig. 1), that is introduced to 
describe the reflections of electrons from the boundary. 

Substituting Eq. (j3Jl into Eq. ((2Jl and taking the initial 
spin polarization in the form P(t = 0) = Pqz, we can 
obtain the following expressions for the components of 
the spin polarization density: 



Px = Pq / G{x -x',y- y', t)^—^ sm{'nr)dy'dx' 

-00 



(5) 



where Go(r,i) = l/(47ri:>i)e-'"'/(4i3t) jg the fuU space 
Green function, D — L'i/{2Tp) is the diffusion coefficient. 



P, = 



(6) 



P. = Po 



G{x ~ x' ,y — y',t)cos{r]r)dy'dx' . (7) 



It is obvious that Py — from symmetry considerations. 
We can not perform the integration in Eqs. © and Q) 
in the general form. However, these integrals can be eval- 
uated explicitly for short times. In this case, since the 
exponents of the Green function effectively cut off contri- 
butions of the terms with large r, the geometric functions 
in Eqs. (|SJ, can be approximated as sin(77r) ~ r/r and 
cos(ryr) w 1 — (?7r)^/2. Then, after the integration, we 
get 



P^ 
Pn 



and 



rjx [ 1 ~ Erf 



— = 1 - rfx^ [ 1 - Erf 
Po ' ' 



j]\j e ""* (8) 



2DtT]^ -1 




V^^dDt / ' 



(9) 



where. Erf [z] is the error function, given by Erf [z] ~ 
(2/ Vtt) / e"* dt. It follows from Eqs. ^ and © that the 



area of nonhomogeneous spin polarization density near 
the edge increases in time and is characterized by the 
length I ~ y/iDt. The behavior of spin relaxation in 
the bulk region is governed by the standard DP theory: 
Px — and the relaxation of Pz is described by the DP 
relaxation time t^ — Tp/ {Lprff' . Let us estimate the lim- 
its of applicability of Eqs. ©, ©. The characteristic 
length scale of the Green functions is \/ADt, while the 
characteristic length of the geometric functions is r]~^ . 
Setting these two lengths equal we find that Eqs. JH)), 
® are vahd while t < 1/{ADt]'^). 

Figure 2 shows the spin polarization density P = 
\J P^ + -P| and its components P^ and P^ {Py — ac- 
cordingly to Eq. ©) as a function of x at several mo- 
ments of time. Spin polarization density represented in 
Fig. 2 was obtained using a numerical integration in Eqs. 
(O and iQ. It is clearly seen the absolute value of the 
spin polarization density P is higher in the edge region. 
It is interesting to note the the higher spin polarization 
density near the edge is largely due to the x component 
of spin polarization density (Fig. 2b). From Eq. Q it 
follows that initial growth of P^ at x — is proportional 
to yi. We also observe formation of a local maximum of 
Pz in the vicinity of the edge. The peak of P^ drifts from 
X — with time. 

In order to check the results of our drift-diffusion ap- 
proximation as well as to find long time behavior of spin 
polarization density, we study the dynamics of spin re- 
laxation near the edge using a Monte Carlo simulation 
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FIG. 3: (Color online) Spin polarization near the edge at 
longer times. Suppression of spin relaxation near the edge, as 
well as spin oscillations in time and space are observed. The 
dotted red line is for the eye. 



scheme that was proposed in Ref. |lOj and subsequently 
used in Refs. |3,[l3- Within the Monte Carlo simulation 
algorithm, the space motion of 2DEG electrons is consid- 
ered to be along classical (linear) trajectories interrupted 
by the bulk scattering events. Our modeling involves 
spin-independent bulk scattering processes, which could 
be caused, e.g., by impurities. For the sake of simplicity, 
the scattering due to such events is assumed to be elastic 
and isotropic, i.e., the magnitude of the electron velocity 
is conserved in the scattering, while the final direction of 
the velocity vector is randomly selected. The time scale 
of the bulk scattering events can then be fully charac- 
terized by a single rate parameter [lOJ, the momentum 
relaxation time, t^. It is connected to the mean free path 
by Lp — vTp. Here v is the mean electron velocity. We 
use reflecting boundary conditions from the edge: the 
longitudinal component of the electron velocity is pre- 
served and the sign of the normal component is changed 
in collisions. The typical number of electrons used in our 
simulations is 2 • 10^. 



We have found that the resuhs of Monte Carlo simu- 
lations at short times are in perfect agreement with pre- 
dictions of drift-diffusion model. Fig. |3| shows spin po- 
larization density obtained via Monte Carlo simulations 
at longer times. At these times the space distribution of 
spin polarization is nontrivial and the slowing down of the 
spin relaxation near the edge becomes more pronounced. 
Fig. indicates the change of the sign of P^ at a; = 0. 
The change of the sign of Px occurs some distance away 
from the edge. Relaxation of spin polarization density in 
the bulk is in agreement with DP theory. 

In conclusion, we have studied the spin relaxation dy- 
namics in a 2DEG near the edge. The obtained result 
indicates that the spin relaxation in the edge region devi- 
ates from the D'yakonov-Perel' theory. Specifically, spin 
relaxation occurs slowly in the edge region and is ac- 
companied by spin oscillations in time and space and 
spin polarization transfer from perpendicular to in-plane 
component. We have derived an analytical formula de- 
scribing initial spin polarization dynamics. Long time 
behaviour of spin polarization density was studied nu- 
merically. Qualitatively, the observed effects are related 



to the lack of translation symmetry in the system. Reflec- 
tions of electron from the boundary partially compensate 
spin precessions due to spin-orbit interaction. Moreover, 
the lack of the actual electron flow from the half-plane 
X < gives rise to non-zero in-plane component of spin 
polarization density near the edge, since this component 
of spin polarization for electrons coming from a; > into 
the edge region is uncompensated. In some sense the 
structure that forms near the edge is similar to long- 
living spin coherence states previously discussed in Ref. 
0. We note that our results are not directly related to 
the spin Hall effect |22|, in which the electron spins be- 
come polarized in an applied electric field. Nevertheless, 
our results could be important in interpretation of spin 
Hall experiments. 
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